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Abstract 

In this paper, we have studied the problem of determining the largest possible set of symmetries for an important example of 
nonlinear dynamical system: the Kuramoto-Sivashinsky (K-S) model in two spatial and one temporal dimensions. By applying 
the classical symmetry method for the K-S model, we have found the classical symmetry operators. Also, the structure of 
the Lie algebra of symmetries is discussed and the optimal system of subalgebras of the equation is constructed. The Lie 
invariants associated to the symmetry generators as well as the corresponding similarity reduced equations are also pointed 
out. By applying the nonclassical symmetry method for the K-S model we concluded that the analyzed model do not admit 
supplementary, nonclassical type, symmetries. Using this procedure, the classical Lie operators only were generated. 
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Similarity Reduced Equations, Nonclassical Symmetries. 



1 Introduction 

As it is well known, commonly the growth process of thin films is from equilibrium and involves the interactions of a 
large number of particles. Notable attentions both in theory and experiment have been devoted to the study of the 
mechanism of thin film growth (see [8] and references therein) . In order to describe the growth process mathematically 
with the partial differential equations of the surface height with respect to time and spatial cooordinates, some models 
have been proposed by researches. These models are based on the collective and statistic behaviour of deposited 
particles. One of them include the Kuramoto-Sivashinsky (K-S) models in which the surface tension, diffusion, 
nonlinear effect and random deposition are considered and so on. 

The K-S equation has been derived e.g. in the context of chemical turbulence [9]. As we declarde above, it is also so 
important because it can describe the flow of a falling fluid film [16]. The solutions of the K-S equation can clarify 
the features about the evolution of surface morphology which is instructive to the investigations of the mechanism 
of thin film growth. 

The two dimensional Kuramoto-Sivashinsky (K-S) model represents a nonlinear dynamical system which is defined 
in a two dimensional space {a;, y}, the dependent variable h = h{x,y,t) satisfies a fourth order partial derivative 
equation which is of the form as follows: 

^ = „V'h-KV^h + X\Vh\\ V = (^,^). (LI) 

where h indicates the height of the interface. The prefactor v is proportional to the surface tension coefficient and 
vV^h is refered to as the surface tension term. The term nV^h is the result of surface diffusion which is because 
of the curvature-induced chemical potential gradient. The prefactor k is referred to as the surface diffusion term 
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and is proportional to the surface diffusion coefficient. The term A| V/i ] represents the existence of overhangs and 
vacancies during the deposition process (when A > 0). The combination of the vV^h and A| V/i | terms models the 
desorption effect of the deposited atoms [8] . 

Some important results concerning the K-S equation are already known. For example, fractal tracer distributions of 
particles advected in the velocity field of the K-S equation was found in [7]. Hong-Ji et al. in [8] have studied the 
evolution of (2-|-l)-dimensional surface morphology in the Kuramoto - Sivashinsky K-S model by using the numerical 
simulation approach. 

The symmetry group method plays a fundamental role in the analysis of differential equations. The theory of Lie 
symmetry groups of differential equations was first developed by Sophus Lie [12] at the end of the nineteenth century, 
which was called classical Lie method. Nowadays, application of Lie transformations group theory for constructing 
the solutions of nonlinear partial differential equations (PDEs) can be regarded as one of the most active fields of 
research in the theory of nonlinear PDEs and applications. 

The fact that symmetry reductions for many PDEs are unobtainable by applying the classical symmetry method, 
motivated the creation of several generalizations of the classical Lie group method for symmetry reductions. The 
nonclassical symmetry method of reduction was devised originally by Bluman and Cole in 1969 [2], to find new exact 
solutions of the heat equation. The description of the method is presented in [6,11]. Many authors have used the 
nonclassical method to solve PDEs. In [5] Clarkson and Mansfield have proposed an algorithm for calculating the 
determining equations associated to the nonclassical method. A new procedure for finding nonclassical symmetries 
has been proposed by Bila and Niesen in [1]. 

Classical and nonclassical symmetries of nonlinear PDEs may be applied to reduce the number of independent 
variables of the PDEs. Particularly, the PDEs can be reduced to ODES. The ODEs may also have symmetries which 
enable us to reduce the order of the equation and we can integrate to obtain exact solutions. In [4], R. Cimpoiasu et 
al. have discussed the classical and nonclassical symmetries of the K-S model in the special case v = k — \ — 1, but 
there are some problems because of some computational mistakes. In this paper, we will try to analyze the problem 
of the symmetries of the K-S equation in a general case by considering all the three prefactors v, k, A. Meanwhile, 
since this paper is the generalization of [4], those computational mistakes exist in [4] will be corrected and some 
further facts are added. 

The structure of the present paper is as follows: In section 2, using the basic Lie symmetry method the most 
general Lie point symmetry group of the K-S equation is determined. In section 3, some results yield from the 
structure of the Lie algebra of symmetries are given. Section 4 is devoted to obtaining the one-parameter subgroups 
and the most general group- invariant solutions of K-S equation. In section 5, we construct the optimal system 
of one-dimensional subalgebras. Lie invariants and similarity reduced equations corresponding to the infinitesimal 
symmetries of equation (1.1) are obtained in section 6. In section 7, is devoted to the nonclassical symmetries of 
the K-S model, symmetries generated when a supplementary condition, the invariance surface condition, is imposed. 
Some concluding remarks are presented at the end of the paper. 



2 Classical symmetries of the K-S Equation 

In this section, the classical Lie symmetry method for the K-S Equation has been performed. First, we recall the 
general procedure for determining symmetries for an arbitrary system of partial differential equations [8,9]. To 
begin, consider a general system of partial differential equation containing q dependent and p independent variables 
as follows 



where u^"' represents all the derivatives of u of all orders from to n. The one-parameter Lie group of transformations 



where = ^[^^o and ip" = ^^\e=a, are given. The action of the Lie group can be recovered from that of its 
infinitesimal generators acting on the space of independent and dependent variables. Hence, we consider the following 
general vector field 



A^(a;,u(")) = 



M = 1, 



(2.2) 



x' = x' + sC{x, u) + 0(£2), u°'=u"+ sip^ix, u) + 0{e^), 



i = l,...,p, a = 1,... 



(2.3) 




(2.4) 
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the characteristic of the vector field V is given by the function 

Q-{x,u^^)) = ^-{x,u)-Y.i\x,u) — , a = l,...,q. (2.5) 



Assume that the symmetry generator associated to (1.1) is given by 

T/ := {x, y, t, h)d^ + f {x, y, t, h)dy + f {x, y, t, h)dt + ^{x, y, t, h)dh (2.6) 
The fourth prolongation of V is the vector field 

= v + if^dh^ + if'dh, + v'dh, + ^^^dh^^ + ^^'dh^, + ... + (2.7) 
with coefficients given by 

= D,Q + C'/ix. + f /ly. + f /it., = D^{DjQ) + i^ha^.j + C'/iy^j + e^/it^j, (2.8) 

where Q = ip — ^^hx — ^^hy — ^^ht is the characteristic of the vector field V given by (2.4) and -D, represents total 
derivative and subscripts of u are derivative with respect to respective coordinates, i and j in the above could be x 
and t coordinates. By theorem (6.5) in [?], the invariance condition for the K-S is given by the relation: 

V(^)[/it + v{h2x + h2y) + Kihix + 2/i(2.)(2„) + h4y) - \{hl + hi)] = (2.9) 
Hence, the invariance condition () is equivalent with the following equation; 

^^^^2. ^ ^2y^ ^ ^y^x ^ 2<p(2-)(2j/) + /j/) _ 2A(<^^/l, + Uy) = 

Substituting () into invariance condition (), we are left with a polynomial equation involving the various derivatives 
of h{x,y,t) whose coefficients are certain derivatives of and (f. Since, and (f depend only on x, 

y, t, h we can equate the individual coefficients to zero, leading to the complete set of determining equations: 





= 0, 




0, 






= 0, 


ey = 


0, 


ilt = 0, = 




= 0, 


- 

S/i — 


0, 


=0, ey= 0, e.t = o 




= 0, 


v/l = 


^0, 


(2A - 2)^, - = 0, (2A - 



By solving this system of PDEs, we find that: 

Theorem 1 The Lie group of point symmetries of K-S (1.1) has a Lie algebra generated by the vector fields V = 

^^{x,y,t, h) = ay + C2t + C3, ^'^{x,y,t, h) = -cax + c^t + ce, 
^^{x,y,t,h) = ci, <f{x,y,t,h) = (^)c2 + (^)c5 + C7. 

and Ci, i = 1, ...,7 are arbitrary constants. 

CoroUciry 2 Infinitesimal generators of every one parameter Lie group of point symmetries of the K-S are: 

d d d d 

Vi = — V9 = — Vq = — Va = — 
ax ay at ah 

Y 9 d ^ ^ d X d ^ yd 
^ ^ dx ^ dy' ^ dx 2Xdh' ^ dy 2Xdh 

The commutator table of symmetry generators of the K-S is given in Table 1, where the entry in the i^^ row and j^^ 
column is defined as [Vi,Vj] = ViVj — VjVi, i,j = 1, 7. □ 
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Table 1 

Commutation relations satisfied by infinitesimal generators 





Vi 


V2 


V3 


V4 




Ve 


Vj 
















V2 


(^)V4 





V2 














Vi 
























Vi 


V2 

























V5 


V2 


-Vi 











Vr 


-Va 









-Vi 





-Vr 








Vr 







-V2 





Ve 









3 The Structure of the Lie algebra of Symmetries 

In this part, we determine the structure of symmetry Lie algebra of the K-S equation. 

has no non-trivial Levi decomposition in the form = t K gi , because Q has not any non-trivial radical, i.e. if r be 
the radical of g, then g = r. 

The Lie algebra g is solvable and non-semisimple. It is solvable, because if g*^^^ =< Vi, [Vi, Vj] >= [g,0], we have: 
gW = [g,0] =< -V^,-l^v^^V^,Vj,-V^ >, and 



5(2) = [0(1), bW] =<^V4>, 



so we have the following chain of ideals g'^^ D g(^) D {0}. Also, g is not semisimple, because its killing form 
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oJ 



is degenerate. 

Taking into account the table of commutators, g has two abelian four and two dimensional subalgebras which are 
spanned by < Vi, 1^2, V3, 14 > and < VQ,Vr >, respectively, such that the first one is an ideal in g. 



4 Reduction of the Equation 

The equation (1.1) can be regarded as a submanifold of the jet space J''(]R^, M). So we can find the most general group 
of invariant solutions of equation (1.1). To obtain the group transformation which is generated by the infinitesimal 
generators Vj = (^jdx + ^fdy + (,fdt + ^Pidh for i = 1, 7, we need to solve the seven systems of first order ordinary 
differential equations. 



dx{s) 
ds 

dyjs) 

ds 
dt{s) 

ds 
du{s) 

ds 



^l{x{s),y{s),t{s),h{s)), 

^fixis),y{s),t{s),h{s)), m- 

^f{x{s),y{s),t{s),h{s)), m-- 

iPi{x{s),y{s),t{s),u{s)), u{0) 



y, 



(4.10) 
(4.11) 



1,...,7. 
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exponentiating the infinitesimal symmetries of equation (1.1), we get the one parameter groups Gk{x) generated by 
Vk for k = 1, 7. 



Theorem 3 The one-parameter groups Gi{t) : M 
table: 



M generated by the Vi, i = 1,...,7, are given in the following 



Gi{s) 


(x, J/, t, u) 


G,{s) 


{x,y,t,u) 


Gs{s) 


(x,y,t,u) 


Gi{s) 


{x,y,t,u) 


G5(s) 


{x,y,t,u) 


Ge{s) 


{x,y,t,u) 


Oris) 


{x,y,t,u) 



■ {x + s. y, t, u), 
' {x,y + s,t,u), 

■ {x,y,t + s,u), 

■ {x,y,t,u + s), 

• (x eos(.s) + y sin(s), y cos(s) — x sin(s), t, u) 

■ {x + St, y, t, ( / ) (ts^ + 2xs + 4uA - 4u)), 

■{x,y + St, t, (^^— ^) (^s^ + 2a;s + 4uA - 4u)). 
where entries give the transformed point exp{tXi){x,y,t,u) = {x,y,t,u). 

These operators which generate invariancc of the evolution equation (2) corre- spond to the following transformations: 
UO generate temporal translation, Ul and U4 Galilean boosts, U2 spatial dilatation and gauge transformation, U3 
a gauge transformation, U5 a spatial rotation, U6 and U7 spatial translations. 

Recall that in general to each one parameter subgroups of the full symmetry group of a system there will correspond 
a family of solutions called invariant solutions. Consequently, we can state the following theorem; 

Theorem 4 If u = f{x,y,t) is a solution of equation (1.1), so are the functions 



Gi{s) .f{x,y,t) = f{x + s,y,t), 
G2{s) .f{x,y,t) = f{x,y + s,t), 
Gsis) - fix, y, t) = f{x, y,t + s), 
G^is) .f{x, y, t) = f{x, y, t) - s, 

Gsis) .f{x,y,t) = /(a;cos(s) +ysin(s),?ycos(.s) - a;sin(s), i), 

2xs + ts^ 

Ge{s) .f{x,y,t) = f{x + st,y,t) - ( )^ 
G^{s) .fix, y, t) = fix, y + st, t) - (^g^) . 



(4.12) 



Thus, from above theorem we conclude that: 

Corollary 5 for the arbitrary combination V = Y^i=i ^ G B) i^e K-S equation has the following solution 

u = e'^"'' fixe"'' +t€5 + ei,?/e^'' cos(e6) + 2:sin(e6) + £2, -yze"'' sin(e6) + 2;cos(e6) + e3,te^"'' + e^) - — (4.13) 

a 

where ej are arbitrary real numbers. 



5 Clcissification of Subalgebras for the K-S equation 

Let A be a system of differential equations with the symmetry Lie group G. Now, G operates on the set of the 
solutions of A denoted by S. Let ,s • G be the orbit of s, and H be an r— dimensional subgroup of G. Hence, H— 
invariant solutions ,s G S* are characterized by equality s ■ S = {s}. If /i £ G is a transformation and s €: S, then 
h ■ is ■ H) = ih ■ s) ■ ihHh~^). Consequently, every iJ— invariant solution s transforms into an hHh~^— invariant 
solution (Proposition 3.6 of [13]). 
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Therefore, different invariant solutions arc found from similar subgroups of G. Thus, classification of iJ— invariant 
solutions is reduced to the problem of classification of subgroups of G, up to similarity. An optimal system of 
r— dimensional subgroups of G is a list of conjugacy inequivalent r— dimensional subgroups of G with the property 
that any other subgroup is conjugate to precisely one subgroup in the list. Similarly, a list of r— dimensional subal- 
gebras forms an optimal system if every 1 — dimensional subalgebra of q is equivalent to a unique member of the list 
under some element of the adjoint representation: ^ = Ad(fif) -1), g £ G. 

Let H and H be connected, r— dimensional Lie subgroups of the Lie group G with corresponding Lie subalgebras f) 
and i) of the Lie algebra q of G. Then H = gHg~^ are conjugate subgroups if and only if ^ = Ad(5') • f) are conjugate 
subalgebras (Proposition 3.7 of [13]). Thus, the problem of finding an optimal system of subgroups is equivalent to 
that of finding an optimal system of subalgebras, and so we concentrate on it. 

5.1 Optimal system of one- dimensional subalgebras of the K-S equation 

There is clearly an infinite number of one-dimensional subalgebras of the K-S Lie algebra, g, each of which may 
be used to construct a special solutions or class of solutions. So, it is impossible to use all the one dimensional 
subalgebras of the K-S to construct invariant solutions. However, a well-known standard procedure [15] allows us 
to classify all the one-dimensional subalgebras into subsets of conjugate subalgebras. This involves constructing the 
adjoint representation group, which introduces a conjiigatc relation in the set of all one-dimensional subalgebras. In 
fact, for one-dimensional subalgebras, the classification problem is essentially the same as the problem of classifying 
the orbits of the adjoint representation. If we take only one representative from each family of equivalent subalgebras, 
an optimal set of subalgebras is created. The corresponding set of invariant solutions is then the minimal list from 
which we can get all other invariant solutions of one-dimensional subalgebras simply via transformations. 

Each Vi, i = 1,...,7, of the basis symmetries generates an adjoint representation (or interior automorphism) 
Ad(exp(£Vi)) defined by the Lie series 

Ad(exp(e.yO.VS) = Vj - s.[V.,V,] + y .[F., [V.,V,]] (5.14) 

where [Vi, V,] is the commutator for the Lie algebra, s is a parameter, and i, j = 1, • • • ,7 ([?],page 199). In table 
2 we give all the adjoint representations of the K-S Lie group, with the (i,j) the entry indicating Ad(exp(eyj))l^-. 
Essentially, these adjoint representations simply permute amongst "similar" one dimensional subalgebras. Hence, 
they are used to identify similar one dimensional subalgebras. 

We can expect to simplify a given arbitrary element, 

V = aiVi + a2V2 + asVs + a^Vi + a^V^ + aeVe + arVr. (5.15) 

of the K-S Lie algebra q. Note that the elements of q can be represented by vectors a = (oi, 07) e M.'^ since each 
of them can be written in the form (6.22) for some constants ai, 07. Hence, the adjoint action can be regarded as 
(in fact is) a group of linear transformations of the vectors (ai, 07). 



Tabic 2 

Adjoint representation generated by the basis symmetries of the K-S Lie algebra 



Ad 


Vi 


V2 


Vs 


Va 


Vi 


Vi 


V2 


Va 


Vi 


V2 


Vi 


V2 


V3 


Vi 


Vs 


Vi 


V2 


Vs 


Vi 


Vi 


Vi 


V2 


Vs 


Vi 


V5 


cos(£) Vi — sin(e)V2 


cos(£)V2 -1- sin(e)Vi 


V^s 


Vi 


Ve 




V2 


Vz+sV, + {^)V, 


Vi 


Vr 


Vi 


V2+(^)V4 


Vs + eV2 + {^)Vi 


Vi 



6 



Ad 


V5 




Vr 


Vi 


V5 + eVi 




Vr 


Vi 


V5 - eVi 


' D 






V5 


Ve - sVi 


V7 - eV2 




V5 


Ve 


Vr 




V5 


cos(e)V6 — sm{e)Vr 


cos(£:)V7 + sin(e)V6 




V5 + SV7 


Ve 


V7 




V5 - eVe 


Ve 


V7 



Therefore, we can state the following theorem: 



Theorem 6 A one- dimensional optimal system of the K-S Lie algebra q is given by 

(1) : V2+aVe, (4): 0^3 + ^7, 

(2) : aV^ + bVr,, (5): V^+aV^ + bV^, 

(3) : aVs + Ve, (6): 0^3 + 14 + 6^5, 



where a, 6, c e R and a^O. 



Proof: F/ : — > g defined hy V 1-^ Ad{exp{siVi)-V) is a linear map, for i = 1, ■ • • , 7. The matrix M? of Ff, 
i = 1, • • • , 7, with respect to basis {Vi, • • • , V7} is 



/I 





















/ 1 





















/ 1 
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1 
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M| = 
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M| = 
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1 










—s 











1 






















1 

















-■< 





1 






















1 







—s 














1 





Vo 

















1/ 


V 








-< 








1/ 


\ 


—s 














1/ 



/I 























-s 



















1 














"A 




i" 


c 














'A 








1 






















1 























1 











Mi = 











1 























1 






















1 























1 
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-s 


Vo 

















1/ 




Vo 
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c J 



/I 
















o\ 


/I 






















1 






















1 





sC 











s 





1 


s^^ 

* 2 
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1 


* 2 




















1 






















1 























1 





s 
















1 


—s 




















1 






















1 





Vo 

















1/ 




Vo 

















1/ 



respectively, where S = sins, C = coss and C = jx- Let V = X^^^i aiVi, then 

F"^ o Fq" o ■ ■ ■ o Fl' : V ^ (5.16) 
(cos(s5)ai - sin(s5)a2 + (C^e cos(s5) - (sr sm{s5))a4)Vi 
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+ (sin(s5)ai +003(55)02 + ((57003(35) + C^e 3in(s5))a4)y2 
+ {sQai + 5702 + as + ^(4 + sl)a4)V3 + 04^4 

+ 1 ( - S2 cos(s5) + Si sin(s5))ai + (s2 sin(s5) + si cos(s5))a2 



+ ((S2 3in(s5) + Si C0S(S5))CS7 + (-S2 €03(55) + sin(s5))CS6) 04 + tl5 - 57^6 + Settrj V5 

+ ^ - S3 cos(s5)ai + S3sin(s5)a2 + (CS3S7 sin(s5) - (ssSq cos(s5) - Csi)a4 + 003(55)06 - sin(s5)a7^ Ve 
+ ^ - S3 3in(s5)ai - 53003(55)02 + ( - (5357003(55) - Cs3S6sin(s5) - Cs2)a4 + sin(55)o6 + 003(55)07^1/7. 

Now, we can 3implify V as follow3: 

If 04 ^ wc can make the coefficients of Vi, V2 and Vq vanish by F^^ , Fg" and F^^ . By setting si = 

*2 = *6 — ^iid 57 = —fj^, respectively. Scahng V if necessary, we can assume that 04 = 1. So, V is 

reduced to the case (6). 

If 04 = and Oi 7^ then we can make the coefficients of V2, V5 and Vq vanish by F^^, F^^ and ^2^. By setting 
55 = — arctan(^), S3 = and 52 = respectively. Scaling V if necessary, we can assume that oi = 1. So, V is 
reduced to the case (5). 



If 04 = Oi = and 02 7^ then we can make the coefficients of 1/3, and V5 vanish by F^'' , F^^ and F^^ . By setting 
57 = — 53 = ^ and 5i = — respectively. Scaling V if necessary, we can assume that 02 = 1. So, V is reduced 
to the case (1). 

If oi = 02 = 04 = and oe 7^ then we can make the coefficients of V7 and V5 vanish by F5 ^ and F^'' . By setting 
S5 = — arctai 
the case (3). 



S5 = — arctan(^) and S7 = fj, respectively. Scaling V if necessary, we can assume that og = 1. So, V is reduced to 



^S6 

6 

respectively. Scaling V if necessary, we can assume that 07 = 1. So, ^ is reduced to the case (4). 



If oi = 02 = 04 = og = and 07 7^ then we can make the coefficients of V5 vanish by Fq^ . By setting 55 = 



If oi = 02 = 04 = 06 = 07 = then V is reduced to the case (2). □ 
5.2 Two-dimensional optimal system 

The next step is constructing the two-dimensional optimal system, i.e., classification of two-dimensional subalgebras 
of 0. This process is performed by selecting one of the vector fields as stated in theorem (6). Let us consider Vi (or 

Vi,i = 2, 3, 4, 5, 6, 7). Corresponding to it, a vector field V = aiVi -\ + 07V7, where Oj's are smooth functions of 

{x, y, z, t) is chosen, so we must have 

[Vx,V]=Wx+wV. (5.17) 

Equation (5.17) leads us to the system 

Cji^ajUk = 'dai + vjai (i = l,---,7). (5.18) 

The solutions of the system (5.18), give one of the two-dimensional generator and the second generator is Vi or, 
Vi,i = 2,3,4,5,6,7 if selected. After the construction of all two-dimensional subalgebras, for every vector fields of 
theorem 6, they need to be simplified by the action of adjoint matrices in the manner analogous to the way of 
one-dimensional optimal system. Thus the two-dimensional optimal system of Q has three classes of g's members 
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combinations such as 

< aiFi + aaVa, + + PsVs >, (5.19) 

< aiVi + a2V3, P1V2 + /32V4 + /33V5 >, 
<Vi+aVz,VT> , <Vi+aV2,V'r > 
<Vi+aiV6,PiV4 + I32V5> ; < Vi + aiVs + a2Ve,Vi > . 

where at, i = 1,2 and /3j, j = 1, ..,3 are real numbers and a is a real nonzero constant. All of these sub-algebras 
are abelian. 

5.3 Three-dimensional optimal system 

This system can bo developed by the method of expansion of two-dimensional optimal system. For this take any 
two-dimensional subalgebras of (5.19), let us consider the first two vector fields of (5.19), and call them Yi and Y2, 
thus, we have a subalgebra with basis {Yi,1^2}, find a vector field Y = aiVi -|- • • • -|- arVr, where aj's are smooth 
functions of {x,y,z,t,u), such the triple {Yi,Y2,Y} generates a basis of a three-dimensional algebra. For that it is 
necessary an suSicient that the vector field Y satisfies the equations 

[Y,,Y] = ^,Y + n7iFi + P1Y2, [Y2,Y] = + ^172^1 + P2Y2, (5.20) 

and following from (5.20), we obtain the system 

q^/J^afc = ^lai + wiP; + piPl, C]kPiak = ^2ai + W2PI + P2PI. (5.21) 

The solutions of system (5.21) is linearly independent of {Yi,!^} and give a three-dimensional subalgebra. This 
process is used for the another two couple vector fields of (5.19). 

By performing the above procedure for all of the two couple vector fields of (5.19), we conclude that Y = l3iYi+(32Y2. 
By a suitable change of the base of q, we can assume that F = 0, so that g is not a 3— dimensional sub-algebra. 
Thus, we infer that: 

Corollary 7 The K-S equation Lie algebra q has no three-dimensional Lie subalgebra. 



6 Similarity Reduction of K-S equation 

The K-S equation (1.1) is expressed in the coordinates [x, y, t, u), so we ought to search for this equation's form in 

specific coordinates in order to reduce it. Those coordinates will be constructed by looking for independent invariants 
{z,w,r) corresponding to the infinitesimal symmetry generator. Hence, by applying the chain rule, the expression 
of the equation in the new coordinate leads to the reduced equation. 

We can now compute the invariants associated with the symmetry operators. They can be obtained by integrating 
the characteristic equations. For example for the operator, H5 := Vg = t-^ — ^x-^ this means: 

dx ^dy ^ (U ^ -2Xdh 

t ~ ~ X ^ ■ ^ 

2 

The corresponding invariants are as follows: z = y, w = t, r ~ h + f^. 

Taking into account the last invariant, we assume a similarity solution of the form: h = f{z,w) — and we 
substitute it into (2) to determine the form of the function f{z,w): We obtain that f{z,w) has to be a solution of 
the following difl'erential equation: 

U + yfzz + kfiz - - ^ = (6-23) 

Having determined the infinitesimals, the similarity solutions zj, wj and hj are listed in Table I. In Table 11 we list 
the reduced form of the K-S equation corresponding to infinitesimal symmetries. 
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Table 3 

Lie Invariants and Similarity Solutions 



J 


Hj 




Wj 




hj 






1 


Vi 


y 


t 


h 


f{z, 


w) 




2 


V2 


X 


t 


h 


f(z, 


w) 




3 


V3 


X 


V 


h 


/(« 


w) 




4 


Vs 


x^ + y^ 


t 


h 


/(« 


w) 




5 


Ve 


y 


t 




/(« 


w) — 




6 


Vr 


X 


t 




/(^, 


w) — 


4At 


7 


V2 + Ve 


t 


y-j 






w) — 


4Ai 


8 


V3 +V5 


x^ + y^ 


t — arctan(|) 


h 


/(^ 


w) 




9 


V3 + V6 


y 


-2x + 


I, , 3tx-t^ 
"■ ^ 6A 




w) — 


6A 


10 


Vs +Vr 


X 


-2y + e 


+ 6A 


/(^ 


w) — 


3ty-t^ 
6A 


11 


Vi + V3 + Vt 


t — X 


y + - tx 


1 x'^ —3x'^ t-\-6xy 
+ 12A 


/(« 


w) — 


12A 


12 


V3+V4+V5 


2 1 2 

X +y 


t — arctan(|) 


h — arctan(-) 




w) + arctan(|) 



Table 4 

Reduced equations corresponding to infinitesimal symmetries 



Similarity Reduced Equations 



U + yfzz + - A/f = 

fw + l^fzz + kfAz - \fz = 

u{f,, + + fe(/4. + 2/..™„ + /4„) - A(/,' + /') = 

+ 4ziy/,, + 4://, + IQkz^U, + GAzf.zz + 32fc/z^ - 4Az/| = 
/» + i^/.. + fe/4. - A/f - = 

/» + I'/z. + fc/4. - A/2 - 2^ = 



2At 

2AaV^ + 2v\{z'^ + z^)U^ + 2k\{z + if Uv: - 2X''{z'' + z^)f^ 
{4k + vz)U^ - Xzfl - AXz^fl + &Akz^f3z + IGkz^U, + {Avz' + 32kz'')U + ^^z^fz + kU^ + z'' U + Skz^f,,^^ = 

16//A/„„ - 4i/A/,, - 64fcA/4» - 32fcA/,,»„ - 4fcA/4z + 16AV^ + 4AVf - w = 

-AuXf,, - IGuXU^ - 4fcA/4. - 32fcA/,,„„ - 64A:A/4.„ + 4AV' + 16A"/^ -w = 

32kX{z^ + z'^)f4w + (96kXz^ + 32fcA)/3» + (IGz^ + 48fcA + 16i/A)/™„ + {96kXz'^ + 32feA)/^z™^ 

+64fcA(a + + 16!/A/,, + (16A - 8Aw)/, - 16X^(z^ + l)fl + {l&vX - \&Xzw)U 

-32X'^zfn,fz + 192kXzf^w^ + 64kXzfzzzw + 32i'Xzfzw + 8vz + 4w^ = 

{4i^z^ + 32kz^)f,, + + 2A«)/^ + 16fe«*/4^ + (^^-^ + 4fc)^,„ + A;/4«, - 4A«Vz + 64^0^3^ 

+4;/zVz + Skz'^f^^^n, - Xzfl -Xz = 



7 Nonclassical Symmetries of K-S Equation 

In this section, we will apply the so called nonclassical symmetry method [2]. Beside the classical symmetries, 
the nonclassical symmetry method can be used to find some other solutions for a system of PDEs and ODEs. 
The nonclassical symmetry method has become the focus of a lot of research and many applications to physically 
important partial differential equations as in [6,11,5,1]. Here, we follow the method used by CAI GUOLIANG et al, for 
obtaining the non-classical symmetries of the Burgers-Fisher equation based on compatibility of evolution equations 
[3]. For the non-classical method, we must add the invariance surface condition to the given equation, and then apply 
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the classical symmetry method. This can also be conveniently written as: 

I^^^^Ai|ai=o,A2=o = 0, (7.24) 
where V is defined in (2.5) and Ai and A2 arc given as: 

Ai := ht + h2x + h2y + /i4x + 2/i(2x)(2s,) + h^y - {hi + hi), A2 := - ^^/i^ - S.'^hy - ^^ht 

Without loss of generality we choose = 1. In this case using A2 we have: 

ht = ^p-i^hx-^''hy. (7.25) 
First, total differentiation Dt of the equation gives 

Dt{ht) = Dt{~i^h2x - i'h2y - khix - '2kh(^2x){2y) - khiy + A/i^ + Xh^y) 

Dt{^ - i^Ux - S.^Uy) = -vh2xt - i^h2yt - khixt - '2khj^2x){2yt) - khiyt + 2Xhxhxt + 2Xhyhyt 

= -UOxxi^P - S.^Ux - i'^Uy) - l^Dyy{(p " ^^Uj; " ^^Uy) " kOxXXxi^fi " ^.^Ux " ^^Uy) 

- kDyyyy{(p " ^^Ux " ^^%) " 2kDxxyy{V " ^^^x " ^^Uj/) + 2 A/l^; Z) ( " ^^Ux " ^^Mj/) 

+ 2XhyDy{(p - (^Ux - S.'^Uy). 

Substituting ^^Uxt and (."^Uyt to both sides, we can get 

^ _^^xx _ ^^yy _ j^^xxxx _ 2k^^^yV _ k^WW + 2A/l^(p^ + 2Xhyipy + S}uxi + S,^Uyt 

Uxxx "I" ^ Uxxy "I" C "^vyx "I" ^ ^yyy + kC '^xxxxx ~\~ k^ Uxxxxy ^" 2k^ Uxxxyy 
~t~ ^ Uxxyyy 2A^ '^x'^xx 2A^ Ux'^xy 2A'Uy^ '^xy 2XUy^ '^yy 

By virtue of 

Dx{ht) = Dx{-vh2x - vh2y - khix - 2kh{2x)(2y) - kh^y + Xh^ + Aft,^) 
Dy{ht) = Dy{-vh2x - i'h2y - kh4x - 2kh(2x)(2y) - khiy + Xh^. + Xh^) 

gives 

hxt ~ ~^hxxx ^h2yx khxxxxx ^ '^kh(2x)(2yx) ^ kh/iyx + 2Xhxhxx + 2Xhyhyx 
hyt ^h2xy ^^yyy kh^xy 2kh^2x}(yyy} ^^yyyyy ~^ 2Xhxhxy ~t~ 2Xhyhyy 

so it gives the governing equation 

ip* + <^2^ + ip^y + (^^^ + 2(^(2r.)(2j,) ^ ^iy _ 2X{ip'=Ux + ^p" Uy) = 

where <^*, are given by 

= Dt{(p - ^^Ux - ^^Uy - ^^Ut) + (^Uxt + ^"^Uyt + ^^Utt = Dt{ip - (^Ux - ^"^Uy) + Uxt + ^^^yt, 
= Dx{ip - (^Ux - ^"^Uy - i^Ut) + i^Uxx + ^'^Uxy + ^^Uxt = Dx{ip - (.^Ux - S.'^Uy) + ^^Uxx + ^"^Uxy, 

= Dyiip - S}Ux - ^'^Uy ~ $^Ut) + ^^Uxy + ^'^Uyy + (^'^ U ly = D I {(f - S}Ux - ^"^Uy) + U xy + Uyy , 

(^^^ = Dxxi^ - ^^Ux - ^"^Uy - ^^Ut) + S}uxxx + S^Uxxy + ^^Uxxt = Dxxi'f " ^^Ux ~ £,^Uy) + S}uxxx + S^Uxxy, 
"P"^ = I^yyif - i^Ux - S.'^Uy - ^^Ut) + £}Uxyy + ^'^Uyyy + Ufyy = Dyyiifi " (^Ux " ^^My) + £}Uxyy + (^Uyyy, 

V ~ I^xxxxi}P C C "^y C ^t) "I" C "^xxxxx "I" ^ '^yxxxx "I" C "^txxxx 

^xxxx{}P C C ^y) ~t~ C '^xxxxx ~^ C ^yxxxx; 

^ yy — D^^yy{{p ^ lia; ^ IVvy ^ t£f ) + ^ 'flxT:ryy ~t~ ^ '^xxyyy ~t~ C '^xxyyt 

— Dxxxxi,^ C ^x ^ ^y) ^~ ^ '^xxxxx ^ ^ ^^yxxxx: 
f^yyyy _ Dyyyy{ip — ^ Wa; — ^ Uy — ^ Ut) + C Uxyyyy + 4 U-yyyyy + C " 



-^yyyyi'P ^ '^x ^ Uy) + ^ "^xyyyy + ^ M 



yyyyy 
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Substituting them into the governing cqiiation, we can get the determining equations for the symmetries of the 
K-S equation. Substituting ^'^ = 1 into the determining equations we obtain the determining equations of the 
nonclassical symmetries of the original equation (1.1). Solving the system obtained by this procedure, the only 
solutions we found were exactly the solution obtained through the classical symmetry approach (theorem 1). This 
means that no supplementary symmetries, of non-classical type, are specific for our model. 

Conclusion 

In this paper by applying the criterion of invariance of the equation under the infinitesimal prolonged infinitesimal 
generators, we find the most general Lie point symmetries group of the K-S equation. Also, we have constructed 
the optimal system of one-dimensional subalgebras of K-S equation. The latter, creats the preliminary classification 
of group invariant solutions. The Lie invariants and similarity reduced equations corresponding to infinitesimal 
symmetries are obtained. By applying the nonclassical symmetry method for the K-S model we concluded that the 
analyzed model do not admit supplementary, nonclassical type symmetries. Using this procedure, the classical Lie 
operators only were generated. 
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